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We consider the Coulomb drag between two metal rings in which the long range Coulomb inter- 
action leads to the formation of a Wigner crystal. The first ring is threaded by an Ahranov Bohm 
flux creating a persistent current Jo. The second ring is brought in close proximity to the second 
and due to the Coulomb interaction between the two rings a drag current .In is produced in the 
second. We investigate this system at zero temperature for perfect rings as well as the effects of 
impurities. We show that the Wigner crystal state can in principle lead to a higher ratio of drag 
current to drive current Jd/Jo than in weakly interacting electron systems. 

I. INTRODUCTION 

Coulomb DragJbfijtsssee n current carrying systems has been studied in numerous papers both experimentally 
and theoreticallyouutllj. The earliest experiments were done in two dimensional systems. In this configura- 
tions the electrons are combined to two-dimensional layers. In one such layer a current is driven and another 
conducting layer is brought to within close proximity of the first. The Coulomb interaction between the 
electrons in these two layers causes a transfer of momentum and the second layer acquires a current. This 
current is referred to as the drag current. 

Recently, the phenomenon of Coulomb drag in one dimcnsionaLsjistems has attracted much attention, 
particularly with regards to nanowires and nanotechnology l3'BEH13ild. The basic description of the one 
dimensional case is identical to that in two dimensions. In all of these systems the time averaged charge 
density of the systems are translationally invariant. The coupling of the two systems is due to charge 
fluctuations in the two systems. 

It is known however that the long range character of the Coulomb interaction can lead to the formation of 
a Wigner crystal in which the electrons become localized and form a periodic lattice. Numerical simulations 
on one dimensional sjjjstems indicate that an arbitrarily weak long range interaction will lead to the formation 
of a Wigner crystalEZEa. This long range nature of the Coulomb force should be apparent for low electron 
densities. 

In this paper we consider the drag between two such Wigner crystals. It is evident that this system has a 
non-uniform charge distributions and we will show that this in principle can lead to a much higher value of 
the drag current. 

II. CLASSICAL DRAG 

It is the aim of this section to analyze from a classical perspective some aspects of current drag that we 
may expect in a quantum mechanical treatment of the Wigner crystal. The Wigner crystal is as a system 
where the electrons are effectively localized, so it is reasonable to expect that some aspects of the current 
drag between Wigner crystals may be exhibited in a classical system. 

Here, we discuss the current drag between two one dimensional wires consisting of classical particles 
interacting via an unscreened Coulomb potential. The two systems are close enough that they are coupled 
by the Coulomb interaction between the wires. This system is illustrated in Figure (nl). To get an idea of 




FIG. 1: Illustration of the current drag setup for considered in this classical analysis. A constant force / drives 
system one. U(x — y) is the interwire interaction and d is the distance of separation. 
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addition, the first chain is driven by a constant force /. In steady state the equations of motion of the two 
chains are given by 

ax\ = — A sin[(27r/a)(xi — £2)] + / 



0LX2 = Asin[(27r/a)(:Ei — X2)] 

Setting <j) = (xi — X2), results in the following equation for the difference between the velocities of the two 
chains. 

2A 

4> = f/a sin(27r0/a) 

a 



with solution 

Tn 



a/2 



l-M sin ( 2 ^/ a ) 



where To is the time for the two chains to slip one complete cycle (distance a/2). 
Carrying out the integral results in 



(00 for / < 2 A/ a (locking) 



This simple analysis shows the interesting result that for suitable conditions of the parameters / and A, the 
two chains will be locked together with the drag current equal to the drive current. 

We now wish to investigate the behavior of this classical system when internal dynamics are included. We 
take therefore as a model of the dynamics of these classical particles the following Fokker-Planck equation 

|p ( f , y, *) = I>(|Jf + j£F p & M + E D (£* + l^fW/) m * t] (1) 

Here, x = {x\,X2...,Xn) labels the positions of N particles in wire 1, y = (t/i, y2---> Vn) labels the positions 
of particles in wire 2, D is the diffusion constant and T is the temperature. P(x, y, t) is then the probability 
of having the N particles in layer 1 at positions x and the N particles in layer 2 at position y at time t. The 
forces present are 

in wire 1 and 

in wire 2. The interaction U is taken to be the unscreened Coulomb interaction between the wires, V is the 
unscreened Coulomb interaction within each wire and / is a constant force driving the particles in wire 1. 
The basic problem is to determine the drag current produced in wire 2 due to the moving charges in wire 1. 

In the numerical analysis of this problem we performed Monte Carlo simulations on discrete systems. 
The conducting systems are partitioned into a lattice with the particles occupying positions at the lattice 
points. The dynamics of the simulation are determined by equation ([!]) which is approximated numerically 
as follows. For each particle, one of two possible random directions is chosen . If the particle lowers its energy 
by moving one lattice site in that direction, then the move is accepted. Otherwise the move is accepted with 
probability e~@ AE , where (3 = l/ksT and AE is the change in the particles energy upon making this move. 
This process is completed iteratively for each particle until a steady state condition is reached. Simulations 
were performed to determine the dependence of the drag current on separation between the two systems. 
The results are shown in Figure (||) 

The plot of Drag Current vs. Separation shows the fact that the sum of the drive and drag currents 
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FIG. 2: Plot of drag current against wire separation for two one-dimensional wires. 



III. WIGNER CRYSTAL RING 



A. Disorder Free Case 



In this section we calculate the non-dissipative drag between two perfect Wigner crystal ri ngs , a t zero 
temperature. The basics of non-dissipative drag effects has been considered in previous papers! 1 ti 15 !. The 
basic system of interest here consists of two metallic rings one of which is threaded by an Ahranov-Bohm 
flux, <£>. This flux shifts the ground state, E g , of the ring to a current carrying state which is calculated 
according to 

T _dE g 

If another ring is brought close to the first one, then the Coulomb interaction will cause a drag current 
to be produced in the second ring. Since this is the ground state of the two ring system this current is 
non-dissipative and persists as long as the magnetic flux is present in the first ring. 

In these works the two rings have uniform charge densities and the drag effect is due to coupling of charge 
fluctuations in each ring. In a system with zero disorder, the relative magnitude of the drag current produced 
is typically a small fraction of the driving current. For two such systems separated- by a distance ~ 200 A 
the drag current is on the order of ~ 10~ 4 times the current in the driving system.li-3 Based on the analysis 
of the last section, we expect that the drag current of two Wigner crystals could in principle be as large as 
the driving current. _ 
I— The model of the one dimensional Wigner crystal that we use here is that proposed in references (E3) and 
(113). In this model the Wigner crystal is viewed as an elastic chain of spinless electrons. In the continuum 
limit, the Wigner crystal of length L, in the presence of an Ahranov Bohm flux, <E>, is described by the 
LagrangiarO 



£ = dx 



171 ( • 2 2c 



2 V) 2 ) 



L$ 



(2) 



where a is the crystal period, m is the electron mass, s is the velocity of sound in the crystal and <&o is the 
flux quantum. The field variable u(x, t) describes the local displacement of the chain at a point x at time t. 
If the field variable is expanded in a fourier series 
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FIG. 3: Schematic illustration of the two ring Wigner crystal setup. The crystal period is a and ui(x,t) denotes 
the local displacement of ring 1 and 112(2;, t) denotes the local displacement of ring 2. Ring 1 has an Arhanov-Bohm 
flux,$i, present which creates a persistent current Jo in ring 1. The interaction between the two crystals, V(ui — U2) 
is a function of the relative displacement of the two cryatals and creates a drag current Jd in ring 2. 



ignored and only the rotation of the ring as a whole is relevant. We wish to study the interaction between 



two such Wigner crystal rings as shown in Figure 
Wigner crystal rings the following Hamiltonian. 



We take as a model for two disorder free interacting 



H 



1 

2M 



-ih 



d_ 



1 

2M 



-ih- 



2n , ^ 

+ V cos — (iti - u 2 ) 
a 



(3) 



where M = Nm is the total mass of each ring, Vb represents the strength of the inter-ring interaction, 
$1 is the flux through ring 1, $2 is the flux through ring 2 and represents the displacement of ring 
1,2 respectively. The first two terms in Equation (||) represent the individual dynamics of rings 1 and 
2 respectively. The third term represents the interaction between the two rings. It is evident that the 
interaction between the chains will be periodic in the relative displacement between the two crystals. For 
simplicity we model this interaction as a cosine term, however the basic features of this analysis should not 
depend on the particular form of the periodic potential used. 



Transforming to center of mass and relative angular coordinates 
have in dimensionless form 



2tt 
L 



(ui — u 2 ) and 8 = f(ui + u 2 ), 



wo 



d 



4 (4 + /3 



q cos(N(f>) 



(4) 



where a = $1 - $2, P = $1 + $2, q = V (MR 2 /h 2 ) and e = E(MR 2 /h 2 ). Our task is now to solve 
equation (0) for the eigenenergy e. 

The hamiltonian, equation (0), separates into two parts H = + Hg with energies e = + eg. The 
ground state wavefunction is a product ^(<fi, 0) = S(8)T(tfi). The equation for the center of mass coordinates 
is readily solved giving 



S{0) 



'2tt 



1 



(0 



{n' + (3) 2 ;n' = 0,±1,±2, 



(5) 



The equation for the relative coordinate wavefunction is 

2 



l g2+2) -<?cos(iV(/>) +e 



T(4>) = 



(6) 
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The energy eigenvalues may be found numerically to any desired accuracy by truncating the set of equations 
for some CK—eaid searching for the value of e which gives a zero value for the appropriate determinant of 
coefficients. c3 

The energy eigenvalues determined in this way, combined with the results of equation (|^) allow the currents 
in each ring to be determined from J\ = and J2 = t^-|$ 2= o-"- The results are shown in Figure (^). 

Figure (Q) shows that the drag current does approach the drive current in magnitude as the interaction Vo 
is increased which is qualitatively similar to the classical result. For large values of the interaction strength, 
q/N > 3 or 

V > 3h 2 N/mR 2 , (7) 

the two Wigner crystals have essentially the same current with each having half the current that one ring 
would have if isolated. 

In order to relate the interaction in equation ([}]) to the true Coulomb interaction we approximate the 
amplitude Vq as follows. We want to calculate the potential felt by an electron a distance d from a finite 
periodic system of length L and period a. We model our system by taking a length L of an infinite line of 
electrons centered at x. This system is shown schematically in Figure (ft). 
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FIG. 5: Schematic model of a periodic system of electrons. The potential on an electron centered at x and a distance 
d from an infinite line of electrons is obtained by taking a finite segment of length L of the infinite chain centered at 
x. 



This potential energy is given by 



n=-±(L/2- x) ^(x-naf+<P 
which can be rewritten as 

r ±(L/2+x) „2 00 
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FIG. 6: Schematic illustration of the two ring Wigner crystal setup. The phase 01,2 = 2zLwi,2(a;, t) where ui,2 is the 
local displacement of the Wigner crystal. The interaction between the two rings is V(fa — fa) =W[1 — cos (fa — fa)]. 
Ring 2 has an impurity of strength V p that couples to the phase fa. Figure (Hb) shows a phase shift in the relative 
coordinate (fa — fa-) of ty relative to Figure (pp.). 



L ' 2 dz' 



-L/2 a \Jz' 2 + d 2 



1=1 



L/2 dz' e 2 {2%l , 

cos z 

- L/ 2 a Vz' 2 + d 2 \ a 



2nl 

COS I X 

a 



The first term is a constant in x therefore we ignore it. For small separation distances d/a < 1, higher I 
terms may be ignored and we keep only the I = 1 term. This yields finally 



V 



, cos (22V) 



L/2 

dz'- , 

Vz' 2 + d 2 



2ir 

COS I X 

a 



We can now identify the amplitude Vq in Equation (^)for a crystal with N electrons as 

4 f L/2 , cos 22V 
Vq = N- / dz' a . 
a Jo y/V^+d 2 



(8) 



4 ,L/ao cos ( 27rz ) 

N- I dr V 

a J y/z 2 + d? 



As mentioned in the introduc tion .for a one dimensional system, a long range interaction will lead to the 
formation of a Wigner crystalEZrL3. The long range nature of the Coulomb force should h& apparent for 
electron densities, n < jq^b where as is the Bohr radius for the particular ring materialo For a typical 
mesoscopic ring radius of 100/im, crystal lattice constant a = lOa^ and as ~ 2A we have L/a > 10 5 so the 
upper limit in the integral in Equation m) may be taken to infinity which yields 



„ 4Ne 2 T ^ .„ d. 

Vo = K 2tt- 

a a 



where Kq is the zeroth order Bessel function. This implies an exponential decrease in the amplitude of 
the interaction with increasing separation between the two Wigner crystals. Comparison with Equation 
shows that the two crystals will be "locked" together at a distance d/a « 4. So two very low density perfect 
Wigner crystals will be essentially locked together when separated by a distance of four times the lattice 
constant or less. 



B. Impurity Effects 



We now consider the effects of impurities on the dynamics of the Wigner crystal. In the presence of an 
impurity, the phase of the Wigner crystal will be pinned. In this case the crystal ring cannot rotate as a 
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Where 

' ' = 8tt 2 1 flt " dx 



ma J C0i 2 2 50i 2 \ 




^ = ^ Nf - si he - W [1 - cos(0 2 )] , (11) 



and the phase 4>\,2 = Ux,i{x, t){2ir/a). 

The first two terms in Equation (J9|) are the Lagrangians describing the dynamics of rings 1 and 2 respefe 
tivcly. The last term describes the interaction between the two rings. We follow the work of Krive et al.EZI 
and introduce the dimensionless parameter 



msa 

which characterizes the magnitude of the quantum fluctuations of the Wigner crystal, and the parameter 

T s = hs 

which characterizes the energy scale in a Wigner crystal of length L. We will be concerned with stiff crystals 
that are weakly fluctuating characterized by a <C 1. r-i 

The basic approach that we will use here is a semiciassicai approximation used by Larkin and Leeta in 
the context of charge density waves and elaborated by KleinertEj for a single particle tunneling through a 
barrier . In this approach the probability amplitude of a soliton tunneling through an impurity is given by 
e A / h where A is the action associated with the tunneling trajectory obtained by minimizing the function 
J dt J dxC{r), where C(t) is the Lagrangian, Equation (O), with t replaced by the imaginary time r = it. 
This approach amounts to approximating the path integral of the propagator by its value at the classical 
trajectory.El 

We now want to answer the following question: given a current in ring 1 that causes a rotation of the 
phase by 2tt, what is the probability that the phase of ring 2 will also rotate by 2ir? It is apparent that there 
are three possiblities: 

i) 0i — > 2-7T, 02 — * 0: the phase in ring 2 will remain pinned and ring 1 will rotate on top of the relative 
potential between the two ending with a relative phase difference (<f>i — <f>2) = 27r. 

ii) cf>i — * 27r, <j>2 — > 2tt: the phase of ring 2 will tunnel through the impurity thus increasing by 2tt and 
ending in phase with <f>i- 

Hi) </>i — > 27r, 4>2 — * —2n: the phase of ring 2 will counter-rotate with that in ring 1 ending with a relative 
phase difference, (</>i — 02 ) = 4-7T. 

Here we have limited our consideration to single tunneling trajectories. In other words, each shift of the 
phase of each ring occurs due to a single tunneling event and ignore all multistep tunneling trajectories. 
Thus we only consider strong coupling between rings, WL/{^j^;) > 1, and strong pinning aV p ^ T s . We 
now calculate the tunneling amplitudes for each case above. 

In case (i), the Lagrangian in Equation (W reduces to 



£=-^ m +* 2 m -wh-cobwo] d2) 




This is simply the problem of a single Wigner crystal ring in the presence of a periodic pinning potential. 
The classical wave equation resulting from this Lagrangian in real time, is 



ma 
8^ 



i-tffa > + Wsin(0i) = 



This is the sine-Gordon equation and it admits the exact soliton solution 

/ x — ut 



0! = 4 tan 1 
where x — —x, t = Lo^t and 



exp ± 
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This solution has the following interpretation: at t = the phase of the crystal is 0, as t increases, the phase 
at x = locally increases to 2n creating a region around x = with phase 2tt. As t increases this region with 
2tt phase propagates symmetrically from x = at speed u towards x = — L/2 and x = —L/2 until the entire 
ring is at <f> = 2tt. The Lagrangian, Equation ( |l2] ) has been considered in the context of charge density wave 
tunneling by Bogachek et alH , and there it is noted that the homogeneous (space independent) soliton 
solution is advantageous in the action against any spatially dependent solution for a finite system. The 
homogenous solution corresponds to a uniform rotation of the Wigner crystal as a whole and has the formEil 



<pi (r) = 4 tan 1 exp(woT). 



(13) 



This solution describes the uniform increase of the phase of the entire ring from to 2tt in a time To ~ l/u>o- 
The action associated with this solution is by Equation (Il2|) 



LU 



WL 



1/2 



(14) 



In case (ii), the relative difference between the two phases does not increase, while the phase of the center 
of mass of the two ring system increases by 2tt. We assume that during the tunneling process the two rings 
remain in phase to minimize the potential energy between them. If we rewrite Equation (^) using (j> — 4>\ — (j>2 
and = 4>i + (j)2, we get the relevant Lagrangian for case (ii) to be 



C = 




1 — cos - 
2 



(15) 



This is the problem of a single Wigner crystal ring tunneling through an impurity. We consider the case 
of strong crystal pinning, aV p T s . 

The case of soliton tunnclmg through an impurity in a strongly pinned Charge Density wave was first 
considered by Larkin and LeeE3. In their analysis, the tunneling process is broken into two stages. The first 
stage consists of a rapid tunneling of a length Iq of the crystal through the impurity in a short time to and 
leaves the crystal in a distorted state of high potential energy. The second stage consists of a slow relaxation 
of the crystal from this distorted state back into a state that minimizes its potential energy. 

To describe the first stage of the process Larkin and Lee proposed the following trial solution to Equa- 
tion @ 



2nt 



for < t < t and \x\ < io- 

It is evident that this solution interpolates between the initial state, 9 = for all x at t = 0, and the 
final state 9 — 2tt for x — at t — 2tt. The second stage of the process takes place away from the impurity 
\x\ > Iq and obeys the equation of motion 



+ s 



2 nil 



0. 



The solution to this stage of the tunneling should obey the periodic boundary conditions of the ring and 
match the [Solution for the first stage at x = —L/2, L/2. This problem has been considered in the work of 
Krive et altZl where the proposed trajectory for the relaxation stage has the form 



7r ± 2 arctan 



(r - t s )s 



The action associated with the tunneling stage of the process is 

h 



A 2 



where are constants. The action associated with the relaxation stage is 



h. ( L \ 
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where C is a constant and C3 is a constant assumed to be of the order unityli3. Equation (16) shows the 



interesting result that the action is only logarithmically dependent on the impurity for strong pinning. 

In case (Hi) above, the center of mass of the two ring system does not move while the phase of the relative 
coordinates increases by Att. The relevant Lagrangian for this system is given by 



C 




d.i 



- V p 6(x) 



1 — cos 



- W 



1 — cos 



(17) 



This is identical to the lagrangian in Equation ( |l5| ) with the addition of the term proportional to W. Since 
the probability for this process goes as e A it is expected that the probability for this case will be reduced 
relative to that in case (ii) by ~ e~ w and can therefore be neglected. 

We are now in a position to analyze the relative magnitude of the drag current Jo to the drive current 
Jo- This ratio is given by 



Jd 
Jo 



en + e~ 



(18) 



It is clear that if \A\\ ^> \A2t\ the two crystal will essentially be locked together. For purposes of 
illustration we consider a stiff crystal with a — .1 and consider values of V p for stiff pinning aV p 3> T s . 

Figure ([?]) shows that the impurity effects are not that drastic in determining the value of W where the 
locking of the two crystals occur. This is due to the fact that the impurity potential enters logarithmically 
into the tunneling action. Comparison with Figure (^b) shows that the distance at which the two Wigner 
crystals become "locked" together is essentially the same as in the impurity free case. The main effect of the 
impurity is then in the magnitude of the persistent current. 

We can estimate the magnitude of these currents by noting that if the two currents are locked together, 
then the magnitude of the persistent current for a given flux $ should be the same as a-sangle Wigner crystal 
ring with twice the mass. This current has been calculated at zero temperature to betll 



Jo ~ Ts 



aVn 



l/c 



e - C3/Q cos I 2tt 



$0 



(19) 



We see that the current in Equation ( |l9|) is much smaller than the disorder free persistent current due to 
the (1/Vp) 1 / 01 factor. 

In summary we have analyzed the Coulomb drag between two one-dimensional Wigner crystal rings. For 
sufficiently large interaction between rings the drag current is essentially equal in magnitude to the drive 
current. For an impurity free ring this "locking" of currents occurs at a separation distance of d ss 4a with 
a the crystal period. We analyzed the effect of a single impurity in the drag ring and found that the two 
currents are "locked" at essentially the same distance and the major effect of the impurity is to significantly 
decrease the magnitude of the persistent current. 
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